Parametrically excited water surface ripples as ensembles of oscillons 
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We show that ripples on the surface of deep water which are driven parametrically by monochro- 
matic vertical vibration represent ensembles of oscillating solitons, or quasi-particles, rather than 
waves. Horizontal mobility of oscillons determines the broadening of spectral lines and transitions 
from chaos to regular patterns. It is found that microscopic additions of proteins to water dra- 
matically affect the oscillon mobility and drive transitions from chaos to order. The shape of the 
oscillons in physical space determines the shape of the frequency spectra of the surface ripple. 

PACS numbers: 47.35.Pq, 47.27Cn, 47.52.+j 



Ripples on the surface of vertically vibrated liquids, 
often referred to as Faraday waves to honor his pioneer- 
ing effort PQ, are commonly used to study generation of 
regular patterns (e.g. [2j [3]) and wave turbulence (e.g. 
[4, 5 ). Surface ripples are usually explained by nonlinear 
interactions of waves in the gravity- capillary range [6j [7] . 

Recent experimental studies of parametric surface 
waves revealed several problems with the wave interpre- 
tation of the disordered surface ripples. Simultaneous 
measurements of the frequency and wave number spec- 
tra [5] reveal poor correspondence between the two: fre- 
quency spectra show distinct multiple harmonics, while 
the wave number spectra are continuous. It has also been 
found that the dispersion relation is satisfied only approx- 
imately. The origin of the universally observed multiple 
frequency harmonics and the mechanism of their spectral 
broadening also remain poorly understood. It has been 
shown that in the time domain, initially continuous waves 
appear to be broken into ensembles of "envelope solitons" 
by modulation instability [8HTU]. What remains unclear 
is how these solitons are realized in physical space. 

In this Letter we show that the surface ripples in a 
parametrically driven system consist of oscillating soli- 
tons, or oscillons, previously found in vertically vibrated 
flows [TTHT5] . We also show that the seemingly chaotic 
state on the water surface can be turned into an ideally 
ordered matrix of oscillons by reducing their mobility. 

We report results obtained in a vertically shaken cir- 
cular container (178 mm diameter, 30 mm deep) with 
monochromatic forcing. The experimental setup and the 
measurement techniques are described in references [8]- 
HD] . Waves are generated in the range of frequencies of 
fo = (10-200) Hz. Figure [TJa) shows a photograph of the 
surface ripple in water, forced at fo = 60 Hz at a peak- 
to-peak acceleration of a = 1.2g. The wave field looks 
random; it is dominated by oscillating "blobs" approxi- 
mately 4 mm in diameter. Blobs move randomly on the 
surface; they merge and collide with other blobs. The fre- 
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quency spectrum of the surface gradient, measured using 
reflection of a very thin (0.5 mm) laser beam off the sur- 
face, is not random. It consists of spectrally broadened 
harmonics, f n = nf /2, Fig.[TJc). 

For the same vertical acceleration of a = 1.2g, the 
surface ripple changes dramatically when a very small 
quantity of a protein suspension (skim milk, gelatine and 
albumin were tested) is added to water, as shown in 
Fig. [ljb) . In this example, we have added 0.5 p. p.m. 
(in weight) of bovine serum albumin to water. Such mi- 
croscopic additions of proteins do not change surface ten- 
sion (a = 73 mN/m) and viscosity (rj = ImPa-s at 20C) 
of water. However, the resulting wave field resembles 
standing waves which form very stable square pattern. 
The frequency spectrum does not change qualitatively, 
but spectral harmonics become very narrow, as seen in 
Fig. 1(c). Red diamonds in Fig. [TJc) show the spectral 
powers contained in each spectrally broadened harmonic 
in water (integral over Af = f n ± 15 Hz). The spec- 
tral energy contained in each harmonic in water does 
not change with the addition of proteins compared with 
the peaks in spectra in water/protein solution. This is 
consistent with the fact that forcing (acceleration) and 
dissipation (viscosity) remain constant. 

The dynamics of the formation of the square pat- 
tern in the albumin water solution is illustrated in Fig- 
ure [2] Shortly after the shaker is turned on, a standing 
wave in the round container forms circular wave fronts, 
Fig. [2|a). Later these concentric circles become mod- 
ulated azimuthally, Fig. |5|b). The standing wave ap- 
pears to be broken into an ensemble of blobs, which 
eventually self-organize into a square matrix in steady 
state, Fig. [5|c). To test if the matrix periodicity is con- 
sistent with the capillarity-gravity wave dispersion rela- 
tion, we perform measurements at different forcing fre- 
quencies. The distances between maxima in the circular 
standing waves agree well with the dispersion, confirm- 
ing their wave nature, Fig. [2|d). However, the distances 
between adjacent peaks in a square matrix noticeably 
deviate from the wave lengths expected from linear dis- 
persion relation uo 2 = gk + (cr/p)k 3 , where k is the wave 
number and p is the density. The discrepancy between 
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the blob periodicity and the dispersion relation suggests 
that the observed pattern is not a standing wave. It is a 
matrix of non-propagating spatially localized oscillating 
surface perturbations. Indeed, the insertion of two metal 
plates into the liquid does not destroy a global pattern 
[Fig. |2je)], as should happen with standing waves sensi- 
tive to the boundary conditions. 

Oscillating blobs, which organize into a square pattern, 
resemble oscillating solitons, or oscillons. Such localized 
oscillatory perturbations of the liquid surface have been 
studied since 1984. The first parametrically driven sta- 
tionary solitonic structures were discovered on the water 
surface in a resonator [11]. Later oscillons, were found 
in granular layers [12 , in thin layers of highly dissipative 
fluids [13], in non- Newtonian fluids [14], and in strongly 
dissipative liquids vibrated at two frequencies [15 . Re- 
cently a non-propagating solitary wave on the water sur- 
face was observed in a very narrow vertically vibrated 
cell below the threshold of parametric excitation by ex- 
ternally perturbing the surface [16]. The dissipation in 
that system was high due to strong wall friction. 

Here we report the first observation of an oscillon in 
a 3D flow driven monochromatically on the surface of a 
deep layer of liquid (larger that the wave length). First, 
we study a single isolated oscillon in a highly viscous wa- 
ter solution of glycerine (77 = 190mPa • s). A shaker is 
driven at the frequency of fo = 30 Hz at the accelera- 
tion just below the threshold of parametric excitation, 
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FIG. 1: (a,b) Photographs of the surface ripples under para- 
metric excitations at 60 Hz and acceleration of a = 1.2g in (a) 
water, and (b) water with a small (0.5 p. p.m. in weight) ad- 
dition of the bovine serum albumin, (c) Frequency spectra of 
the surface gradient measured using reflection of a thin laser 
beam off the fluid surface in water (red) and with albumin 
added (black). 




FIG. 2: Flow development during the startup phase in the 
albumin solution at (a) t = 1.5 s, (b) t = 2.5 s, and (c) 
t = 6.0 s. (d) Spatial period of a square pattern as a function 
of driving frequency: solid squares correspond to the distance 
between the maxima along the diagonal of the matrix; open 
squares correspond to the distances in rows/columns. Red 
circles show spatial periods of the concentric circles at the 
early stage of the ripple evolution. Solid line shows dispersion 
relation of the capillary- gravity waves, (e) A square pattern 
is not substantially affected by the insertion of metal plates 
into the established matrix. 

a = 3.1g. The liquid surface stays flat until it is per- 
turbed externally (e.g. by a falling droplet), as seen in 
the frame sequence of Fig. |3ja,b) obtained from the fast 
video (300 fps) movie [17 . A localised oscillating struc- 
ture is sustained by the forcing for a very long time (over 
several minutes). 

The shape of the oscillon was studied by illuminating 
its cross section by a thin green laser sheet in combina- 
tion with a fluorescent dye (Rhodamine B) which was 
added to the liquid. The laser sheet produces a sharp 
image, as shown in Fig. |3|c). The profile of the surface 
contour is well approximated by the hyperbolic secant fit, 
Fig.^d). This is in agreement with the oscillon shape 
in other systems (e.g. [H]). The frequency spectrum 
of the surface gradient of a single subcritical oscillon is 
shown in Fig. j3^e) for the excitation frequency of /o = 60 
Hz. Remarkably, this spectrum shows the subharmonic 
of the excitation frequency as well as its harmonics. Mul- 
tiple harmonics are observed universally in all experi- 
ments with parametrically excited waves. Their presence 
is usually attributed to wave- wave interactions or even to 
the energy cascade in wave turbulence (e.g. |5]). Here we 
show that multiple harmonics is a feature of an individual 
oscillon. 
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FIG. 3: (a,b) A single oscillon excited at the subcritical ac- 
celeration (a = 3.1g) in viscous water solution of glycerine by 
external perturbation at 30 Hz. (c) A snapshot of the oscil- 
lon profile illuminated by a vertical laser sheet after adding 
fluorescent dye to the glycerine-water solution, (d) surface el- 
evation of the oscillon is well approximated by the hyperbolic 
secant fit. (e) Frequency spectrum of the surface gradient of a 
single subcritical oscillon. (f) Surface ripple on the glycerine- 
water solution above critical excitation (a = 3.3g). (g) Fre- 
quency spectrum of the surface gradient of the supercritical 
ripple. The excitation frequency is 60 Hz. 

When the acceleration is increased above the threshold 
of the parametric excitation, multiple oscillons are gen- 
erated. They form a square matrix, Fig. [3jf), identical 
to those observed in water with added protein (Fig. [2]). 
Their frequency spectra measured on the surface of the 
oscillon matrix, Fig.[3|g), are very similar to those of in- 
dividual oscillons, Fig. |4jc). Thus, frequency spectra of 
the surface motion are determined by individual oscillons 
and not by collective effects. 

At higher forcing, oscillons become more mobile in hor- 



izontal plane. This increased mobility gradually leads 
to a less regular pattern, as seen in a fast video [18]. 
A difference in the mobility of oscillons determines the 
spectral widths of the frequency harmonics. Since the 
shape of an oscillon in physical space is given by the 
hyperbolic secant, h(x) ~ sech(ax), its slow (compared 
with the period of oscillation) movement about the ob- 
servation point should lead to a modulation envelope 
h(t) ~ sech[ir / (at)] in the time domain. This, in turn, 
should lead to a spectral line shape in a frequency spec- 
trum given by the squared hyperbolic secant function 
of frequency, E(f) ~ sech 2 [b(f — / )] [8 j. The analy- 
sis of the surface ripples of Figs, [I] and [3] confirms these 
expectations. Fig. [4^ a) shows shapes of the spectral 
lines (first subharmonic) in pure water and in the albu- 
min solution (same conditions as in Fig. [I]). A higher 
mobility of oscillons in water leads to the line shape 
E(f) ~ sech 2 [b(f — fo)]. The addition of albumin, re- 
duces the line width to the instrumental limit determined 
by the measurement time. In the time domain, the higher 
oscillon mobility in water leads to the observation of the 
"envelope solitons" previously reported in [HI [10]. The 
results presented here reveal that the hyperbolic secant 
envelopes result from the shape of the individual oscil- 
lons in physical space. The sech analytic fit to the time- 
domain signal is illustrated in Fig. Qb). In contrast to 
this and consistent with the narrow spectrum, stationary 
solitons within the matrix [Figs. [TJb), [2^c) , [3^f )] gen- 
erate non- modulated harmonic signals seen in Fig. [4jc). 
The results similar to those obtained in the water-protein 
suspension have been found in a viscous glycerine-water 
solution, as illustrated in Fig. [4] This suggests that the 
same effect, namely the reduction of the oscillon mobility 
(possibly due to the change in their short-range interac- 
tion) is achieved in Newtonian (glycerine) and possibly 
non-Newtonian (water-protein) liquids. 

Oscillons result from the instability of the surface 
waves, as illustrated in Figs. [2|a-c). This instability in 
the capillary wave range is modulation instability [10], a 
counterpart of the well known Benjamin- Feir instability 
of gravity waves [19 . In its nonlinear stage, this insta- 
bility leads to the breaking of waves into ensembles of 
oscillons, Fig. [2] A similar process is observed in the 
gravity wave range in large wave tanks, where the cross 
wave instability destroys planar wave fronts [20]. The 
connection between cross waves and the oscillon forma- 
tion has been suggested earlier [21]. Since oscillons are 
formed in the process of breaking of the surface wave, 
they inherit their size from the mother wave such that 
their diameter and the period of patterns that they form 
roughly follow the dispersion relation, Fig. [2^d). 

Individual, or subcritical, oscillons shown in Fig.[3]were 
observed in a viscous glycerine water solution. They have 
also been observed in water with added protein suspen- 
sions, though their amplitudes were lower there. Stronger 
external perturbations in water lead to the onset of multi- 
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FIG. 4: Shapes of the spectral lines on the water surface 
illustrated in Fig. [I] (a) Frequency spectral line of the first 
subharmonic (30 Hz) in water (triangles), water-glycerine (di- 
amonds), and in water-albumin solutions (squares). Solid line 
shows sech 2 fit. Waveforms of the surface gradient at the fixed 
point on the surfaces of (b) water, and (c) water-glycerine and 
water- albumin solutions. 

pie oscillons on the surface and to a pattern formation. A 
mechanism through which microscopic additions of pro- 
teins (typically less than 1 p. p.m.) dramatically change 
the mobility of oscillons is not understood. It is well 
known, however, that small additions of high molecular 
weight polymers have striking effects on the stability and 
pattern selection in turbulence, for example in Couette- 
Taylor flow [22] . This is attributed to different viscoelas- 
tic properties of these liquids. As we show here, oscillons 
are observed in both Newtonian (glycerine) and in weakly 
non- Newtonian fluids. The effect of proteins on the pat- 
terns of surface ripples should be further investigated. 

Our results suggest that patterns and ripples are made 
of oscillons which interact on a short range. There are 
two types of subharmonic oscillons which differ by a 
phase shift of l//o- As in other oscillon systems (e.g. 
[12]), the like-phase oscillons show repulsive interaction, 
while oscillons of the opposite phase attract and bind. 
Such interaction forces lead to a variety of effects, in- 
cluding the formation of nontrivial patterns, e.g. square 
matrices in a circular container, Figs. [2jc) , |3jf), as well 
as stripes (when the matrix symmetry breaks along rows 
or columns). 

It should be noted that the secant hyperbolic solitons 
are known as a limiting case of cnoidal waves. A transi- 



tion from a periodic cnoidal wave to a sech solitary struc- 
ture determined by the elliptical parameter k {k < 1 for 
waves and k = 1 for solitons) has been discussed in rela- 
tion to the oscillon theory [23] . 

The interpretation of the spectral broadening as the re- 
sult of the increased mobility of oscillons may be relevant 
in other wave systems as well. Exponential frequency 
spectra are observed universally in many parametrically 
excited systems, such as the magnetic spin waves excited 
in ferrites [24 , nonlinear waves in optical fibres [25] , Ra- 
man fibre lasers [26] and in magnetically confined plasma 
[27] . Modeling of the parametrically driven damped non- 
linear Schrodinger equation [28] shows the development 
of nonlinear localized structures, such as oscillons in dis- 
sipative nonlinear systems. 

Finally, in his seminal paper pQ, M. Faraday mostly 
refers to " crispations" made of "heaps" rather than 
waves. 

This work was supported by the Australian Re- 
search Council's Discovery Projects funding scheme 
(DP110101525). 



[1] M. Faraday, Philos. Trans. Royal Soc. of London 121, 
319 (1831). 

[2] M.C. Cross and P.C. Hohenberg, Rev. Mod. Phys. 65, 
851 (1993). 

[3] A. Kudrolli and J. P. Gollub, Physica D 97, 133 (1996). 
[4] W.B. Wright, R. Budakian and S.J. Putterman, Phys. 

Rev. Lett. 76, 4528 (1996). 
[5] D. Snouck, M.-T. Westra and W. van de Water, Phys. 

Fluids 21, 025102 (2009). 
[6] J. Miles and D. Henderson, Annu. Rev. Fluid Mech. 22, 

143 (1990). 

[7] M. Perlin and W.W. Schultz, Annu. Rev. Fluid Mech. 

32, 241 (2000). 
[8] H. Punzmann, M.G. Shats and H. Xia, Phys. Rev. Lett. 

103, 064502 (2009). 

[9] M. Shats, H. Punzmann and H. Xia, Phys. Rev. Lett. 

104, 104503 (2010). 

[10] H. Xia, M. Shats and H. Punzmann, FPL 91, 14002 
(2010). 

[11] J. Wu, R. Keolian and I. Rudnick, Phys. Rev. Lett. 52, 
1421 (1984). 

[12] P.B. Umbanhowar, F. Melo and H.L. Swinney, Nature 

382, 793 (1986). 
[13] O. Lioubashevski, H. Arbell and J. Fineberg, Phys. Rev. 

Lett. 76, 3959 (1996). 
[14] O. Lioubashevski, Y. Hamiel, A. Agnon, Z. Reches and 

J. Fineberg, Phys. Rev. Lett. 83, 3190 (1999). 
[15] H. Arbell and J. Fineberg, Phys. Rev. Lett. 81, 4384 

(1998). 

[16] J. Rajchenbach, A. Leroux and D. Clamond, Phys. Rev. 

Lett. 107, 024502 (2011). 

[17] http: //people .physics . anu.edu.au/~mgsll2/movies/ 
| 30Hz_sub_7_07_2011.mp4| 

[18] http: //people .physics . anu.edu.au/~mgsll2/movies/ 
ripple_glycerine . avi 



5 



[19] T.B. Benjamin and J.E. Feir, J. Fluid Mech. 27, 417 
(1967). 

[20] S. Lichter and L. Shemer, Phys. Fluids 29, 3971 (1986). 
[21] A. Larraza and S. Putterman, J. Fluid Mech. 148, 443 
(1984). 

[22] A. Groisman and V. Steinberg, Phys. Rev. Lett. 77, 1480 
(1996). 

[23] J.W. Miles, J. Fluid Mech. 148, 451 (1984). 

[24] I.V. Krutsenko, V.S. L'vov and G.A. Melkov, Sov. Phys. 



JETP 48, 561 (1978). 
[25] J.N. Kutz, C. Lynga and B.J. Eggleton, Opt. Express 13, 
3989 (2005). 

[26] S.A. Babin et a/., Phys. Rev. A 77, 033803 (2008). 

[27] D.C. Pace, M. Shi, J.E. Maggs, G.J. Morales and T.A. 

Carter, Phys. Rev. Lett. 101, 085001 (2008). 
[28] X. Wang, Physica A 154, 337 (2001). 



